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We examine the prospects for LHC discovery of SU(2)l singlet vector-like quarks that obey Min- 
imal Flavour Violation (MFV) and are consistent with lower energy phenomenology. We study 
models where the vector-like quarks have the same quantum numbers as ur or da, allowing mixing, 
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which generally leads to significant low energy constraints. We find that at leading order in the MFV 

OX) 

expansion there are two naturally phenomenologically viable MFV models of this type when the 
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Weyl spinor components of the vector-like quarks are flavour triplets. We examine direct production 
bounds, flavour and electroweak precision data constraints for these models and determine the cross 
section for allowed values of the model parameters at LHC. For the models we identify as naturally 
phenomenologically viable, large amounts of parameter space afford a significant early discovery 
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reach at LHC while being consistent with lower energy phenomenology 
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In addition to the known standard model (SM) quarks, models beyond the SM generally include new 

in 

matter content. Such extended matter content is strongly constrained by the requirement of gauge anomaly 
cancellation JT). In recent years, vector-like quarks (for reviews see G10) have been studied since the 
. £h possibility of vector-like SU(2)l singlet quarks is a particularly simple and minimal anomaly free extension 

X 

of the SM. From a top-down model building perspective, vector-like matter of this form is also interesting. 
Extensions of the SM involving new vector-like matter have been explored in the context of SUSY to 
address the little hierarchy problem [4—7], and triplets of vector-like quarks can also appear in grand unified 
models and string compactifications that embed the SM into an Eq group [8]. 

The current phenomenological constraints on vector-like matter come from direct production bounds 
at the Tevatron and flavour physics. Tevatron studies exclude a heavy t' with SM like couplings at 95% 
CL up to 335 GeV H and a heavy b' with SM like couplings at 95% CL up to 338 GeV [10;]. Flavour 
constraints are also significant for vector-like quarks. The mixing of n vector-like quarks with the SM 
quarks causes the 3x3 SM CKM matrix to become non-unitary as it is extended to a (3 + n) x (3 + n) 
unitary matrix. Non-unitarity of this form is tightly constrained. By mixing with the vector-like quarks, 



effective tree level flavour changing neutral currents (FCNC) could also be induced for the SM quarks (see 
01 for recent studies). Although flavour constraints are stringent due to the lack of any clearly statistically 
significant evidence of non-SM FCNC's (see, e.g., ifTTl for a survey), these constraints do not provide a 
direct mass bound, as the mixing factors are generally unknown and can be chosen to be small. 

Thus, it is interesting to examine vector-like quark models within the framework of Minimal Flavour 
Violation (MFV) |[T2l - fl6l where new physics effects are flavour conserving apart from the effects of flavour 
breaking by insertions of Yukawa matrices, as in the SM. MFV is also predictive in terms of the allowed 
representations and couplings of new physics. For example, it has been shown that only a single flavour 
singlet scalar representation is allowed by MFV besides the Higgs IfTTl (see lfT8l - l211 for recent studies). 
When the scalar is allowed to transform under the flavour symmetry more representations are allowed 
[22], but MFV remains a predictive formalism that constrains the allowed mass spectra and couplings. 
Similarly, for the vector quark models we study, the strength of the vector-quark quark mixings are no 
longer uncorrelated free parameters in the theory due to MFV. 

In this paper we examine vector-like quark models whose Weyl spinors transform as flavour triplets when 
MFV is an approximate global symmetry. We will show that at leading order in the MFV expansion only two 
models exist that are naturally phenomenologically viable; one with a triplet of charge Q = —1/3 vector- 
like quarks and one with a triplet of charge Q = 2/3 vector-like quarks. By naturally phenomenologically 
viable we mean that the new flavour changing effects in the model are suppressed, in agreement with 
experiment, not simply through the choice of parameters but in a way that follows from the group structure 
and representation content of the theory. This flavour naturalness criteria is a generalization of the Glashow- 
Weinberg criteria for natural FCNC suppression |23l and is not guaranteed by MFV alone. The continuing 
success of the SM in the flavour physics program, with ever more flavour changing decays and flavour 
oscillation measurements probing the weak scale and consistent with the SM, strongly implies that models 
which satisfy this naturalness criteria are promising models to search for in the LHC era. 

The outline of this paper is as follows. In Section II we present our procedure for solving the mass 
spectrum of the models and determine the two models that are naturally viable. We then determine the 
allowed parameter space considering relevant collider, flavour and electroweak precision data (EWPD) 
constraints in Section III. We then turn to the expected cross sections at LHC for the allowed remaining 
parameter space of these models in Section IV Finally, we comment on similar flavour triplet MFV models 
of this form when my > TeV. 



II. MFV VECTOR-LIKE QUARK MODELS 

The SM has a global U(3) 5 flavour symmetry that is only broken by the SM Yukawa interactions. 
Decomposing the abelian quark subgroup of this symmetry as 

G F = SU(3) U ® SU(3) d ® SU(3) Q , (1) 

one can formally restore this subgroup of the flavour symmetry by treating the Yukawa matrices as spurions 
that transform as 

<W ~ (1,3, 3), U ~(3,1,3), (2) 

such that MFV forbids FCNC's at tree level. The flavour symmetry is broken by insertions of the spurions 
9u9u or g u g u when SU(3)tj r or SU(3)q l indices are contracted. We use notation where each insertion 
is accompanied by (the unknown constant) rfr and choose r]i ~ r\ <C 1. Except for the top, all Yukawas 
are already small so the symmetry breaking effects are highly suppressed. For the top, our perturbative 
expansion in r\ is appropriate under this assumption. The naturalness of this assumption depends on the UV 
completion of the models under study and specifying a UV completion is beyond the scope of this work. 
For approaches to MFV without this assumption see irT4Hl6l . In our analysis, due to this assumption, the 
effects of mixing are the dominant constraint on the mass scale of the model. 

In this paper, we are considering models whose left- and right-handed (Weyl spinor) components trans- 
form as flavour triplets. Similar MFV models have been examined before in both the quark and lepton sector 
|[24l|25l . Both of these works assume that the bare mass parameters are far larger than the electroweak scale 
v, which can be appropriate for > TeV scale vector-like quark masses. In this paper, we directly solve 
the models at leading order in the MFV expansion without this assumption and study the constraints and 
prospects for these models in much greater detail. This turns out to be instructive as we argue it singles out 
two models that are naturally phenomenologically viable. We will discuss our approach in some detail for 
an extension of the SM with vector-like down quarks. The discussion for vector-like up quarks is similar. 

A. Vector-like down quark models 

The Lagrangian including the right- and left-handed chiral fields of the vector-like quarks is the follow- 
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c d = Q L iipQ L + d R i]pd R + v£iipvg + v£iipv£ 



+ 



4mf V[ d R + 4™ d 2 v£v£ + 4 (^ ) Q L HV% + g d Q L Hd R + h.c. 



(3) 



The transformations allowed under Gf are listed in Table I for the new couplings and fields. Note that the 
couplings K,f are matrices whose transformation law depends on the model and is such that the Lagrangian 
is invariant under Gp. MFV requires that these matrices be Yukawas or identity matrices. The factors mf 
set the scale of the couplings. The fields V&, V£ are triplets under SU(3) C , singlets under SU(2)l and have 
hypercharge —1/3. As in the SM, when MFV is imposed, Ql transforms as (1, 1, 3) while <1r transforms 
as (1, 3, 1) under Gf- 





SU(3) Ur x SU(3) Dr x SU(3) Ql 


model 


4 


K d 


K d 
K 3 


V* 


Vr 


I 





(1,1,1) 


(1,3,3) 


(1,3,1) 


(1,3,1) 


II 


(1,3,3) 


(1,1,1) 


(1,1,1) 


(1,1,3) 


(1,1,3) 


HI 





(1,3,3) 


(1,3,3) 


(1,1,3) 


(1,3,1) 


IV 


(1,1,1)* 


(1,3,3) 


(1,1,1) 


(1,3,1) 


(1,1,3) 


V 


(1,3,3) 


(3,1,3) 


(3,1,3) 


(1,1,3) 


(3,1,1) 



TABLE I: Representations of Vf, Vr allowed by MFV. We have used the fact that in models I, III the Vj^ has the same 
transformation as <Ir under the flavour and gauge groups to rotate to a basis in V® — Dr space where the mixing 



term vanishes [24]. We also note that model IV predicts non-hierarchical corrections to SM masses due to mixing. 
These corrections require extra tuning to cancel against the SM masses when mf, mf 3> v. This was indicated with 
a star appended to the relevant coupling. 



B. Solving vector-like down quark models 



We now diagonalize the Lagrangians corresponding to the models listed in Table I. There are two forms 
of mixing now to contend with when vector-like quarks are added to the SM; mixing that in the SM corre- 
sponds to the transformation between the weak eigenstates and the mass eigenstates, and mixing between 
the vector-like quarks and the SM quarks. We diagonalize in a two-stage procedure. The fields that corre- 
spond to the initial Lagrangian terms in Eq. ([3]) are designated as the unprimed fields. We introduce unitary 
transformations to primed fields (which would be mass eigenstates in the SM if no vector-like quarks were 
present) and we introduce analogous unitary transformations for the vector- like quarks: 

d R = U(d, R) d' R , d L = U(d, L) d' L , 
u R = K(u,R)u' R , u L =U(u,L)d' L , 
V$ = U(V, R) V$ , V? = U(V, L) Vf . (4) 



We choose them in such a way that the initial Yukawa matrices are diagonalized by taking 1 

U\d,L)g d U{d,R) = * , U\u,L)g u U{u,R) = ^^ . 

v v 

We write the Lagrangian in terms of the primed fields and the mass matrix M d 

T 



(5) 



Cr 



vi> 



M d 



d' R 



(6) 



where 



M d 



M d 



(7) 



rr4tf{d,L)i4U{V,R) 
mjU*(V, L) KJU(d, R) m^U^V, L) n$U{V, R) 

We designate the true mass eigenstates of the SM quarks and the vector-like quarks as double-primed 
fields. The primed basis is now purely unphysical for the vector-like quarks; it does not correspond to the 
eigenstates of any interaction. Thus we are free to choose U(V, R), U(V, L) in order to diagonalize each 
3x3 block of the Hermitian matrices in M d as any rotation of the vector-like quark field basis is allowed. 
This can be accomplished (up to insertions of Vckm) by choosing the U(V, R/L) to be the appropriate 
U(u/d, R/L) or 1 as MFV forces the nf to be (undiagonalized) Yukawa matrices or the identity matrix. 

We determine the M d M d and M d M d matrices in terms of the SM parameters in M. d , At„, and Vckm- 
The 6x6 matrix M d can always be diagonalized by introducing appropriate rotations of the primed fields 
to the mass eigenstate basis: 



d'L 

yd,, 



u 



d' 



Vf 



d'k 

v£» 



w 



d' R 

vi> 



(8) 



where U, W are the unitary transformations diagonalizing M d M d , M d M d , respectively, and D = UM d W^ 
is diagonal. Writing these transformations as 



U 




W 



W21 U>22 



(9) 



where Uij,Wij are 3 x 3 diagonal matrices, the general form for the kinetic terms is 
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3 J/i - 



vk 



-U- 



21" 



92 I gl sin 



e_w\^ 



(J2 sin e w i 
3 ^ 



V? 



+ f 2 u L V C KM^ + [(«n - l)d" L + u 21 V}!] + f 2 [(n n - l)d'i + u 2X V'{\ V CKM ty u" L 

(10) 



u u u 21 ^+A{d!' L $V>l + V£U'L) , 



Although Mi is a real diagonal matrix and v is the vacuum expectation value of the real component of the Higgs doublet, 
until the further mixing effects between the vector-like quarks and the SM quarks are taken into account one cannot identify the 
elements in this matrix in terms of the physical quark masses as yet. 




which will be used to derive the constraints for the vector-like quark models. The up quarks do not require 
any further rotation so that u' L , R = u" L , R . 

Note that the SM CKM matrix is defined as Vckm = U^(u, L) U(d, L) and 

= U{u,L) I Uh \- (ID 

V VcKM<i' L J 

The CKM matrix in these models corresponds to the difference in rotating the up-type quarks in the SU(2)l 
doublet compared to the down-type quarks (which are an admixture of d" L and V£) from the unprimed weak 
basis to the double -primed mass basis. The resulting 3x3 SM CKM matrix is now non-unitary due to the 
correction of the ul W + d^ + h.c. coupling that depends on (tin — 1). An advantage of our approach to 
diagonalization is that the resulting non-unitary 3x3 SM CKM matrix has a clear decoupling limit. 

Using this procedure to solve the Table I models (and a similar procedure for the vector-like up quark 
models) singles out one model in each case that is more naturally phenomenologically viable (at leading 
order in the MFV expansion). We now discuss in some detail this vector-like down quark model before 
turning to the effects of flavour breaking. 

1. A viable vector-like down quark model 

In model I we have nf = 0, k\ = 1, nf = ga- Without loss of generality we have chosen U(V, L) = 
U(V,R) = U{d, R) to work in a basis where the 3x3 subblocks of Eq. dTl) are diagonal. We find the 
diagonal form of the squared mass matrix 




D 2 = | v u/ | , (12) 



?7ln 



where C{ = y2 -^- — ^ and we neglected terms O(ej). From the normalized eigenvector matrices one 
identifies the U and W transformations as 

l ■ J 1_t ~ e \ 2 | , W-={ I -^ I . (1.3) 

u 




' 2 / \ ^3 

where we neglected terms O(ef). Note that £j = mf/v. We do not assume that v -C mf 23 ora large 
hierarchy of masses mf S> m\. We are interested in solutions where the elements of the diagonalized M® d 
(given by m®) are such that m? <C v thus e» <C 1. We note that no further source of CP violation is present 
due to the final rotation that diagonalizes the masses, and the SM masses are identified identically with m® 
in this model. We also note that rotations of V R L can eliminate explicit CP violating phases from the new 
Lagrangian terms in Eq. Q in model I, and that here the indices on the £» are treated as labels coincident 
with the flavour index, not true flavour indices, i.e. they are not contracted. 



This model has a number of features that distinguish it from the remaining models of Table I and make 
it naturally phenomenologically viable: 



(i) The non-unitarity of the CKM matrix is due to a correction of the form 

(14) 
This leads to small deviations from the SM CKM; the largest deviations are in Vtb, V u b, V c b, which 



VcKM = V CKM 
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are less precisely measured, and are proportional to m 2 /v 2 . As we discuss in Section 
deviation from unitarity, being contributed to only the diagonal terms of the square of the CKM 
matrix, removes contributions to meson mixing observables. 

(ii) The modification of the SM couplings of the down quarks due to mixing with the vector-like quarks 
is proportional to 

WO^g)^. (15) 

MFV does not guarantee small mixing effects on the SM quark couplings to the W^ , Z. This model, 
and a similar vector-like up quark model, do have small mixing effects without the mass scale of the 
vector-like quarks being my » v. 

The final phenomenologically relevant couplings are those involving the Higgs that are generated from 
terms initially involving vector-like quarks. Up to order 0(e;), 

C V h = e i a(d'L) l h(Vi"y + e^2(Vi' / ) t h(d'[)\ (16) 

where £2 = — • The SM Yukawa coupling receives a small correction 0(e?) which we have neglected. 
We have also neglected the effects of spurion breaking of the U(3) 5 flavour symmetry but have explicitly 
included the flavour indices. We now determine the flavour breaking effects up to leading order in 77 and the 
dominant decay widths of this model. 

2. Flavour breaking and decays in model I 

The dominant flavour breaking effects come from insertions of the g\ g u or g u g\ spurions. The kinetic 
terms have suppressed flavour indices that one can contract with spurion insertions. Canonically normaliz- 
ing the fields after such insertions ensures that only interaction terms in the Lagrangian that are not bilinear 
in Ql , ur or du will receive corrections of this form. Thus, the only terms that receive flavour breaking 
corrections are 

5C = 6 4 Mr,) 6 (d'Di h (v£"y + e % hirj) 6 (v£")i h (d'ly . (i?) 



Taking into account the canonical rescaling of the Ql field one finds /i (77) = 1 + 77 ^f £^3. These flavour 
symmetry breaking terms will not be the focus of our analysis since they do not provide clearly dominant 
constraints on the mass scale of the theory (they depend on the unknown symmetry breaking parameter 77 
and we assumed 77 <C 1). Our focus will be on the larger flavour diagonal mixing effects with the SM quarks 
that lead to significant mass constraints independent of 77. 

The Vji dominantly decay through Higgs interactions while the V£ dominantly decay through charged 
current interactions to SM up quarks which can afford correlated collider signatures. The largest decay 
widths are given by 
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d\2 



gh 2 j\(ycKM 
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in 



d\3 
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64 7T 



e]e 2 fM 2 

32 7T 
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(m$) 2 



1 + 2 



in 
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(m't) 2 



in. 
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d\1 



(18) 
(19) 
(20) 



C. Vector-like up quarks 



The Lagrangian including the fields VR, V£ is 

C u = Q L i$Q L + u R i$u R + VEi$VE + V%ilt>V R l 



+ 



K u x m\Viu R + K u 2 m u 2 VZV R + «3 



% - !m ^ Ql^VI + guQi^uR + h.c. 



(21) 



The fields VR, V£ are triplets under SU(3) C , singlets under SU(2)l and have hypercharge +2/3. The 
allowed representations are listed in Table II. 





SU(3) Uh x SU(3) Da x SU(3) Ql 


model 


k\ 


K u 


K u 
K 3 


vi 


VR 


VI 





(1,1,1) 


(3,1,3) 


(3,1,1) 


(3,1,1) 


VII 


(3,1,3) 


(1,1,1) 


(1,1,1) 


(1,1,3) 


(1,1,3) 


VIII 





(3,1,3) 


(3,1,3) 


(1,1,3) 


(3,1,1) 


IX 


(1,1,1)* 


(3,1,3) 


(1,1,1) 


(3,1,1) 


(1,1,3) 


X 


(3,1,3) 


(1,3,3) 


(1,3,3) 


(1,1,3) 


(1,3,1) 



TABLE II: Representations of V£, VR. Model IX predicts non-hierarchical corrections to SM quark masses, this was 
again indicated with a star. 



D. Vector-like up quark Lagrangian construction 

If we assume a mild hierarchy mi > ITI3 such that a is still a good expansion parameter for the top 
quark, the vector-like up quark models are nearly identical to the corresponding down-type models except 
for modifications 2 in the kinetic Lagrangian corresponding to having Q = + 1 and T^ L = + \ . The kinetic 
part of the Lagrangian is 



Vffi+ffl -// /* ,// 



C km = C^ + K 1 -l)^^^< + ^' i$-W?&t + 



2gi sin 9w rf i 2 g2 sin 9w 



3 ^ ' 3 



V" 



+ n 



L 



i0 + ( u 2 21 ^M _ 2jejtote\ g + 2 92 sme W j i 



V'r' 



+ f 2 [(un - lK + u 21 V£] V CK M^ + d!' L + f 2 d" L V& KM ty [(«n - 1)< + UaiVl'] 
+ n llU2 i^P K^' + y£%£) . (22) 

As before, we must allow for insertions of both g u g u and g u g u anywhere we have contractions of SU(3) UR 
and SU(3)q l indices, respectively. Thus, insertions of f(rj) occur more frequently in the up-type models 
and provide non-trivial differences between the phenomenology derived in the down-type models. 

1. A viable MFV vector-like up quark model 

The viable MFV up quark model is model VI. The analysis of this model proceeds as in Section II B 
1 with the appropriate replacements. The only difference in the analysis comes from the flavour breaking 
effects. Similarly to model I, we have the following Higgs interaction terms 

C Vh = u 6 / 2 (r?) (ul)i h {VtT + et 6 f 2 (v) (V#")i h «)* . (23) 

2 

The flavour breaking corrections are given by f 2(1]) = 1 + 3 77 ^f ^3 and the partial widths are obtained 
directly from Eqs. (18)-(20) with the appropriate replacements. 

E. Remaining vector-like flavour triplet models 

We have emphasized that the models which we study in detail are naturally phenomenologically viable 
under the assumption 77 <§C 1- We have also solved the remaining models in Tables I, II and summarize 
their mass spectrum and mixing angles in the Appendix. These other models become more viable as the 
mass scale my increases and mixing effects decrease in the decoupling limit, along with LHC discovery 



One replaces e.g., d R — > u R , U(d, L) — > U(u, L), Vckm — > Vqkm' !m< ^ -^2 ~~ > -M" 



10 

potential. One can always choose the parameters in Eqs. Q, ( [2J"| ) to reduce mixing effects for the models 
in the Appendix, violating our naturalness criteria. Such a tuning is also not protected by any symmetry 
against radiative corrections. In the remainder of the paper, we will focus our attention on the two naturally 
phenomenologically viable models. 

III. CONSTRAINTS ON VECTOR-LIKE QUARK MODELS 

We will study EWPD constraints, flavour physics constraints and current collider constraints on mod- 
els I, VI. We begin with the current collider constraints. In each case we analyze the relevant constraint 
assuming that the SM is extended by model I or model VI. 

A. Collider constraints 

1. Vector-like down quarks 

For model I we use the model-independent results from ifTUll constraining deviations from SM predic- 
tions of pp -> Z + (> 3j) at the Tevatron. The addition of singlet vector- like quarks does not supply a 
positive contribution to the T EWPD parameter that can raise the fitted Higgs mass value. As the decays 
h —> ZZ are highly suppressed for the best fit Higgs mass values (recent studies find Mh = 96l 24 GeV 
[26]) only the pair production of Vl will yield significant branchings to Z + (> 3j). Due to this we 
rescale the simulated production cross section (which has appropriate acceptance and phase space cuts for 
the analysis ifTOlO shown in Fig. 38 of this study by \. 

We also rescale to correct for the branching ratio to the constrained final states. The experimental study 
assumed BR(6' — > Z + (> 3j)) = 1. The effective branching ratio for the produced V£ to produce the 
signal (P cS ) is 

/3 cff ~ BR(Vg ->Z<P L )x BR(y[ -> W + u l L ) x BR(W + -> hadrons) 

+ BR(Vl -> ZaP L )x BR(Vl -> W~ v} L ) x BR(W~ -> hadrons) 

+ BR(V[ -»• Zd j L )x BR(VI -)• Z d{) x BR(Z -»■ hadrons) . (24) 

From the PDG we determine BR(Z /W^ — > hadrons) and we assume the produced hadrons initiate jets 
that will pass the triggers of [10]. From Eqs. (18-20) we also determine 



br(v[ -, z d {) * - 2 ( l + J )m L > ^> 
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where we have neglected corrections of 0{m^/rriy), 0(m 2 z /my). Only retaining contributions that are 
not CKM or light quark Yukawa suppressed we find /3gg ~ 0.38 using couplings and masses defined at 
/i = mz- Rescaling the LO a(pp — > b' V) curve of ifTOl in Fig. 38 by -^/3 e ff> we obtain a 95% CL lower 
bound for the vector-like down mass of 200 GeV for model I. 

Recent studies [27] searching for b' b' — > (tW T ) (tW^) improve this bound, assuming the decay prod- 
ucts pass the corresponding triggers. To use this study the effective branching ratio is given by 

f3 cS ~ BR(Vl -> W~ t L ) x BR(y L 3 -». W + i L ) ~ 0.44 . (27) 

Rescaling the results in Table II of IT271 by ^/3 e g we obtain a 95% CL lower bound of 260 GeV. 

2. Vector-like up quarks 

For model VI we use the results of [9] searching for pair produced t' i! that decay into excess iJftx+jets 
events from i! — > Wq decays at the Tevatron. This study selects for one and only one isolated muon or 
electron with Ej- or Pt respectively greater than 25 GeV as a trigger lepton. The remaining trigger requires 
$T > 20 GeV and at least four jets with Et > 20 GeV and \r/\ < 2.0. Again we rescale a(t' i') given in 
[9] by 1/2 as only the left-handed vector-like quarks decay into final states that contain W^. 

We also rescale to correct for the appropriate branching ratio in model VI to the triggered-on final states. 
In this model the effective branching ratio for one of the produced V£ to give the signal (fit) is given by 

ft ~ BR(y[ -»■ W + dl) x [bR(VJ? -»■ W~ di) x BR(T^" -»■ hadrons) 
+ BR(Vjj -»■ W + di) x BR(V^ ^ Z u{) x BR(Z -> hadrons 
+ BR(t^' -> W" (? L ) x [BR(V^' ^ Zu{)x BR(Z -^ hadrons) (28) 

~ 0.61 . 

Rescaling the NLO a(pp —> i! i') curve in Fig. 2 of (9J by -f-(3t we have conservatively obtained a 95% 
CL lower bound for the vector-like up mass of 325 GeV for model VI from the upper limit of the expected 
95% CL region. 

B. Electroweak precision data constraints 

1. EWPDfit 

EWPD constraints are weak as the singlet vector-like quarks do not directly break custodial symmetry 
and for models I, VI contributions to the fl^ with more than one mass scale in the loop are suppressed by 
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€i. In considering EWPD we neglect NP effects that are suppressed by e, and neglect the effects of spurion 
insertions that break the flavour symmetry proportional to rj. For model VI we assume a mild hierarchy 
itt-2 > TUs so that the top Yukawa e, can also be neglected. With these assumptions the contribution to the 
vacuum polarization of the Z (proportional to g^ v ) from a singlet vector-like quark of charge Q is 

n ^ (g2) = gh™ 2 °wN c N f Qig^ ^^^ _ ^ + q2 _ .^ + q2)Bo{q2 ^ ml ^ y)] (29) 

This result is given in terms of the Passarino-Veltman functions ll28ll with standard definitions. All other 
vacuum polarizations can be determined in terms of this result through taking appropriate derivatives and 
rescaling the couplings. We then construct the EWPD parameters 12914311 We use STUVWX ll32ll33Tl as 
the constrained masses are expected to be in the 100 GeV range. The results for models I, VI are reported 
in Table III. 



STUVWX fit results 


model 


68 % CL 


95 % CL 


model 


68 % CL 


95 % CL 


I 


82 GeV 


81 GeV 


VI 


147 GeV 


102 GeV 



TABLE III: The results of an STUVWX fit as defined in [20| where the CL regions are defined for six parameters 
through the cumulative distribution function. The constraints on the vector-like up quarks are stronger as Q = 2/3. 



2. Rh constraints 



The oblique EWPD fit [20] used in the previous section did not include deviations in the parameter 
Rt, = T{Z — > bb)/T(Z — > hadrons). Writing the coupling of the Z boson to the quarks as 



£; 



g\ + g\ Z* q^(f L , q P L + fn, q Pr) q 



(30) 



(3D 



where Pl/r = (1 =F 75)/2, the tree level couplings to the Z are given by 

fl q = T*- S m 2 e w Q, fl q = - S m 2 6 w Q. 

Deviation in the left-handed coupling, due to mixing with the vector-like quarks, is given by fi t f, = 7^6 + 
8 fi^b ■ This leads to a deviation in R^ that can be approximated by ||34l 



<Ji2b~2i#(l-i#) 



Uh? + U°R, b ) 2 



-0.78df L , b 



(32) 
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For the vector-like down quarks the largest anomalous contribution is given by a tree level shift in the 
Zbb coupling. In model I, a positive deviation SfL,b = e|/2 is predicted. As the measured value of Rf, at 
the Z pole and the SM predictions Il35l are 

R rneas _ q.21629 ± 0.0066 , #f M = 0.21578 ± 0.00010 , (33) 

the la bound on the deviation of Rf, from its experimental measurement due to the SM and this deviation 
6fL,b is given by —1.5 x 10~ 3 < 5fh,b < 2 x 10~ 4 . This translates into the following bound on the 
parameters in model I 

For the vector-like up quarks, determining the contributions to Rf, involves modifications of the SM 
diagrams for the top contributions to Rf, and new diagrams where Z — > t V^ in the loop. The required 
calculation has been performed in sufficient generality before and we use the results of [36] for the large 
m 2 , my 3> m^ limit. We retain only the effects due to (V£) 3 that are proportional to m 2 /v 2 . Recalling 
that for model VI the third generation singlet quark has a mass rriy = (m^) 2 (1 + e|), the contribution to 
Rb is given by 

2 



16 it sin 2 6w \ m 2 v 



3{m 2 v -m 2 ) + 2 m 2 v m 2 t /m 2 w I ' m v \ _ m 2 

1 1 o ( 2 J 2 

vfiy-ml V m t ) m w 



.(35) 



C. Flavour constraints 

MFV forbids FCNC contributions to meson mixing at leading order in the MFV expansion that have 
been discussed in the literature J3]|37|. However, the effects of SM quarks mixing with virtual vector- like 
quarks and anomalous SM couplings can still contribute to deviations from the SM phenomenology of 
meson mixing measurements. We first consider these constraints before going to b — > s 7 constraints. 

1. Vector-like down quarks: Meson mixing 

Consider the additional contribution to the effective Hamiltonian for \SS\ = 2, K° — K° mixing, which 
can be characterized by an extra contribution to the dominant operator 

^|<55|=2 = (Csm + C v ) (d L Y s L ) (d L j u s L ) + h.c. . (36) 
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As the virtual quarks are charge +2/3, the potentially largest flavour changing effects for K° — K° 
mixing in model I (when r\ <C 1 and top Yukawa breaking of the flavour symmetry is neglected) come 
about through the non-unitarity of the effective CKM matrix. The anomalous Wilson coefficient (following 
the notation of El) is given by 

Gp a 

nix/ 

a=c,t 



C v 



8 J2 X a B (x a )X sd + X, 



y/2 47r svnOw 
where the effective Inami-Lim function [38, 39] and the non-unitarity of Vckm is given by 



(37) 



B (x a ) = - 



+ 1 \y? log(x a ) 



X. ***>ql IXq -L 






and we have defined x a = {m^) 2 /m^, A" d = V* s V a d- Note that the anomalous Wilson coefficient is SM 
CKM and Yukawa suppressed as expected due to MFV. We observe that in model I non-unitarity effects 
of this form on K°, B®, B° mixing vanish as the deviation from unitarity is proportional to off-diagonal 
elements of the (decoupled) square of the SM CKM matrix J2 a=u c t A" d . This effect is easily seen in a 
simple example with a unitary matrix P and a non-unitary matrix Q where 

( A b\ ( X \ 

P=\ , Q =\ \ , (39) 



A B \ 


, Q = 


/ X 


C D ) 




\ Y 



then 

f 1 + 2X + X 2 \ 

(P + PQ)t(P + PQ)= . (40) 

y 1 + 2Y + Y 2 J 

The deviation from unitarity in the square of the effective CKM matrix is only in the diagonal entries and 
does not effect K°, B®, B mixing through terms proportional to X s d- This form of non-unitarity has been 
termed pseudo-orthogonality in some past studies of vector-like down quarks HOI from Eq compactifica- 
tions. We adopt this nomenclature for this property of the effective CKM matrix. 

Measurements of D° oscillations, however, proceed through charge —1/3 quarks and are sensitive to 
virtual contributions from the singlet quarks in model I. The Lagrangian term of interest is 

C = 92 ( m ^) 4 {V C km){ W + (Vl)j + h.c. , (41) 

\ v m-2 / 

where the index on rrij is again a label coincident with the flavour index and is not contracted. The above 
CKM entries correspond to the SM entries as the vector-like quarks decouple in the limit m^ — > oo. These 
CKM entries are equal to the measured CKM values up to consistently neglected higher order (m'j) 2 jv 2 
terms in model I. Considering CKM and Yukawa suppression, the vector-like quark contribution is domi- 
nated by the contributions from the third generation vector-like quarks and the corresponding contribution 
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to xd (we use the definition of this mixing parameter given in ||4T|) is 

v 2 G 2 F m 2 w f 2 D M D 2 ^( m bmi\ A 
x v D =* % oV B D [V cb V uh ) ri(m c ,m w ) f f{x v ) , (42) 

3 7T Z 1 D \ V TTlo / 



where xy = ray/m^, and f(xy) — ► xy [1 + 6 log(xy)] in the large xy hmit, which is applicable consid- 



ering collider constraints discussed in Section [III A| We use the values of the CLEO-c determination 11421 of 
f D = 222.6±16.7+2i MeV and the lattice calculation of B D = 0.83 reported in R3l The renormalization 



group running for a LO calculation is given by B4l 

6/25 



ri (,,M) = (^Y 2/7 (^) 6/2 7^ 

\a s (m t ) ) \a s (m b ) ) \ a s [p) 



(43) 



The remaining parameters we take from the reported results in the PDG. The measured value of xd from 
the Belle collaboration in the analysis of D° — > K$ vr + tt~ |45l is given by 

x D = (0.80 ±0.29 ±0.17) x W 2 , (44) 

while fits to the HFAG database (without this measurement as a prior) give xd = 8.4^ 3 ' 4 x 10~ 3 [41] from 
the short distance SM OPE contribution. Long distance contributions have been estimated to be of the order 
of 10 -3 ll46l but are difficult to reliably calculate. The contribution of the vector-like down quarks to this 
quantity is far too small to be detected due to Yukawa and CKM suppression. 

2. Vector-like up quarks: Meson mixing 

The effect of charge Q = 2/3 vector-like up quarks on D° mixing vanishes due to the pseudo- 
orthogonality nature of the modified CKM matrix. However, these vector-like quarks can contribute to 
K°, B°, B® meson mixing through virtual contributions in box diagrams. We ignore rj corrections to the 
Higgs coupling in these diagrams as they represent small Yukawa suppressed symmetry breaking terms. We 
only retain the contributions from (V^) 3 whose effects are proportional to m 2 /v 2 . Using the results of ||3) 
the new contributions to M^ for K° — K° mixing are given by 

Mf 2 = — ^ W 1n K o [ X V W So(xy) +2\ c \y 7] c y S (x c , Xy) + 2 X t Xy T] tV S (x t , Xy)] . 

12 7T 



The functions Sq(x), Sq(x, y) are the usual Inami-Lim functions given by 

Ax — llx 2 ±x 3 3x 3 

4(1 -x) 2 2(1 -x 



S ^ X ) ^ —^ =yj 9M_^3 l0 S X 



3xy xy{x 2 -8x + 4) xy(y 2 -8y + A) 

S ° {X ' y) = "4(x -!)(!/-!) + Hx-l) 2 (x-y) l ° gX + 4(y - l) 2 (y - x) l ° gV " (45) 

The contribution of this expression to Aex is 

|AeA-| = C e B K e 3 Im [e 3 X 2 rjy S (x v ) + 2 A c X t rj c y S (x c ,x v ) + 2 X 2 rj t y S (x t , xy)] . (46) 
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Recall £3 = — "'J" 3 . We take NLO values for the approximated QCD running from the top mass which 
gives rjv = 0.58, r\ c y ~ r\ ct ~ 0.47 and rjtv ~ Vtt — 0.57 from O |47j. Using the measured values 

m K = 497.6 MeV, f K = (156.1 ± 0.8) MeV, (AM K ) exp = (3.483 ± 0.006) x 10~ 12 MeV one obtains 

3.65 x 10 4 . 



C f 



G 2 F F 2 K m K M^ 



(47) 



6y/2n 2 AM K 

Further, lattice QCD ll48l gives the input Bk{2 GeV) = 0.54±0.05. Considering that current measurements 
find \ex\exp = (2.229 ± 0.010) x 10~ 3 while recent CKM fits lead to the SM prediction \tK\theory = 
(1.8 ± 0.5) x 10~ 3 , we find the la bound - 0.07 x 10~ 3 < |A e K \ < 0.93 x 10~ 3 . 

We have also determined the constraints in model VI from B° — B° mixing and B® — B® mixing from 
the results of [3]. We find that these constraints are negligible when compared to kaon mixing and the other 
constraints discussed. 



3. b — > s 7 constraints 



Recall that the effective Hamiltonian (neglecting the light quark masses) in the SM is given by 

2 G F 



n c $ 



V2 
2G f 



r errib 



8=1 



V tb V t *A(x t ) 



y/2 "»'»"^' Ll6vr 2 
where the Inami-Lim function for the operator O-j above is 



s L a^b R F llv + 



(48) 



A{x) 



x 
3 



+ 



1 



+ 



2 (x-1) 2 (x-1) 2 (x-1) 



+ 



+ 



9 



+ 



(x-1) 2 (x - l) 2 (x - 1) 



+ - 



x 



3 

log x — - 



logx 



(49) 



(x-1) 4 b 2 (x-1) 4 
and Xj = m 2 /m 2 v . For new contributions to b — > sj in linear MFV, the largest effect of the vector-like 
down quarks of model I is the modification of the SM W^ coupling, while for vector-like up quarks the 
largest modification comes about through the third generation virtual vector-like quark in the usual loop 
contribution to Oj. In both models the largest contribution comes from the third generation vector-like 
quark. 

For model I the largest new contribution to O-j is 



C 



v 



V2 



V tb V t *A(x t ) 



«jx 2 

rribm^ 



while for model VI the largest new contribution is 



C 



v 



±Gj_ 
V2 



V tb V*A(x v ) 



vrrin 



vrrio 



(50) 



(51) 
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One can use the results of [50] for the contribution of such a Wilson coefficient to BR(S 

X s 7)Sv>l.6GeV given by 



BR(B -> X s 7 )^>i.6GeV = 3.15 ± 0.23 - 8.0 



y/2C% d 



4G f V tb V t 



ts 



Comparing to the current world experimental average fl5T 



BR(B -► X s 7) By >i.6GeV = 3.55 ± 0.24+°^ ± 0.03 , 



/2 C u,d 

we obtain a 1<t bound of —0.17 < ( An ,/,„ 



< 0.07. 



(52) 



(53) 



D. Combined constraints 



The constraints on the model parameters mi and mj, that we have derived are presented in Fig. 1. 



Constraints on model I 



Constraints on model VI 




400 600 

m 3 [GeV] 



20 30 

m 3 [GeV] 



FIG. 1: The consttaints on model I (left) and model VI (right). The shaded regions are excluded by the various 
labeled observables considering a \a deviation in the labeled measurement adding the SM and the vector-like quark 
contribution; for EWPD a 68% CL exclusion region is shown, while for Tevatron constraints the 95 % CL exclusion 
region is shown. 



Shown are the 68% CL regions from the EWPD fit, the 95% CL constraints from direct searches at the 
Tevatron utilizing the b' V — > (tW T ) (tW^ decays for model I and £$t+ jets constraints from if — > Wq 
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decays for model VI. Also shown is the la bound for Kq — Kq mixing, R b and b — V sj. Note that the 
parameter space shown for model VI is characterized by m-i > m^ which is consistent with the assumed 
mild hierarchy allowing the vector-like quark Lagrangian construction without retaining all orders in €3. 
Fig. 2 shows how the flavour and R\, constraints are relaxed for models I and VI when larger deviations are 
allowed by comparing the la and 2a regions. 



Rb and flavour constraints on model I 



Rb and flavour constraints on model VI 




400 600 

m 3 [GcV] 



> 

P, 600 







/ A(lcr) 






/ 




/ 


,'' Ka-Ko (Id) ; 




1 / 






,-' v ~----__ 


s^ K„ - K„ (2a) 






-^rf^ 


i^pj^ 


i?t(2cr) 



20 30 

m 3 [GeV] 



FIG. 2: Comparison of the la (dashed line) and 2a (solid line) meson mixing and Rb constraints for models I (left) 
and model VI (right). 



For comparison, to illustrate the relative flavour unnaturalness of the remaining models consider models 
III and VIII, whose mixing angles and mass spectrum are given in the Appendix. (We compare these models 
simply for ease of comparison as only two Lagrangian parameters are present for these models.) Models 
III and VIII do not satisfy our flavour naturalness criteria. One has to consider Lagrangian parameters 
yml + mjj ^ 10 6 GeV for both models so that the lightest vector-like quark is > 100 GeV considering 
direct production bounds. The vector-like quark masses would then be in a pattern mf : m^ ■ mf given by 
~ 10 2 : 10 3 : 10 5 GeV in model III; in model VIII the pattern would be similar with the heaviest vector-like 
quark being ~ 10 6 GeV. Only the phenomenology of the lightest vector-like quark states would be readily 
accessible at LHC or the Tevatron. 

It is convenient to change to a polar coordinate system when determining the constraints on this model, 
where r ~ ym|~+~m| and m?, = iri2 tan</>. Viable parameter space requires r > 10 6 GeV and 0^1, 
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where m l v ~ t/mj r. Expanding in large r and small <fi one finds 



2r 2 +n 2 ^ 1 _ f r +" 



Using this parametrization the constraint from Rf, for model III is given by 



r 2 + v 2 



< 1.5 x 1(T 3 , (55) 



2r 2 

which gives < 0.055, while the direct production bound at the Tevatron is given by r > 7.3 x 10 6 GeV. 
The Lagrangian parameters must be chosen consistently with these strong constraints for this model to be 
viable and any LHC vector quark signal in early runs would have to correspond to this small subset of the 
parameter space. 

For model VIII the correction to fa & is 



a / mtf — m% 



5f L ,b = ~ r^~ v 2 r +3 log 

16 7r sur 0jy V m w 



my 

2 

mf 



a (f> (3 mf — 4m v m 2 + m v ) a 4> [(m v — Sm^) m 2 + 3m v m 2 ^ 
16 7T sin 2 9w (tti 2 / — m 2 ) m^ 8 it sin 2 9\y ("Jy — w- 2 ) "z 2 ^ 



7T4" 1 
2 



(56) 



The vector-like up quark collider constraint from Vl —> W qi decays [9] directly applies for the lightest 
Q = 2/3 vector-like quark giving my > 325 GeV. This translates into a constraint r > 2 x 10 7 GeV. The 
allowed values of (j) are <fi < 2 x 10 -2 for my = 325 GeV which monotonically decreases to <p < 2 x 10~ 3 
for my = 1200 GeV. Again, the Lagrangian parameters must be chosen consistently with these strong 
constraints for this model to be viable. 

IV. PRODUCTION AT LHC 

The discovery of vector-like quarks at LHC has been studied extensively in the literature 
Generally, the discovery signatures rely on the decay of the vector-like quarks producing a W or Z boson 
giving leptonic tags as well as SM quarks that initiate jets. Studies at LHC are likely to follow the Tevatron 
studies of this form (9l [TO] |27l but also have the opportunity to apply new theoretical approaches such as 
employing jet mass Il57ll58ll 

The results of these various studies are applicable to the MFV models studied in this paper with the ap- 
propriate rescaling of the decay widths in terms of the known mixing parameters. However, as most of these 
studies have employed simulations of LO QCD, in this section we present the production cross sections for 
the allowed remaining parameter space in Fig. 1 using the analytic inclusive NLO QCD production results 
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of 11591 which include the partonic production channels 



q + q -»■ v L/R v L/R + x , 
9 + 9^ v L / R v L/R + x , 
9 + q^ Vl/r Vl/r + x . 



(57) 



We correct this result with a factor of Nj /2 because of the three flavours of the vector- like quarks and the 



chiral factor, as only Vl will typically give gauge bosons that can be triggered-on (see Section III A for 
details). The effective cross section is 

N f 



O-cff = -y (o-qq + <J gg + Ogq) 



(58) 



We show the sum of these inclusive cross sections in Fig. 3 (left) determined using the NNLO PDF's of 
MSTW ||60| and four-loop running /3 function results of [61] to run from the reference value a s (Mz) = 
0.1135 ll62l . We vary the renormalization scale of the evaluation between my jl < fi < 2 my to define a 
scale dependent error on the production cross section. The results are shown in Fig. 3 (left) and the effective 
a is presented in Table IV for -y/s = 7, 10, 14 TeV. We additionally show the various contributions of the 
partonic production channels to the inclusive pair production effective cross section for ^fs = 10 TeV also 
in Fig. 3 (right). 




m v [GcV] 



FIG. 3: Left figure: The effective NLO QCD pair production cross sections for y/s — 7, 10, 14 TeV, where the error 
band is defined through the variation of the renormalization scale my/2 < /i < 2 my. Right figure: The various 
contributions to the effective NLO QCD pair production cross sections for ^/s = 10 TeV. The solid line is the sum of 
all contributions. The short dashed line is inclusive pair production through g g —> VlVl, while the long dashed line 
is inclusive pair production through qq — > Vl Vl, which is highly suppressed due to the antiquark PDF suppression. 
Production through qg — >• VlVl dominates for high masses. 



21 





NLO 
°eff 


my [GeV] 


^=7TcV 


V^ = lOTcV 


V^ = 14TcV 


332 


10i| Pb 


27±l 2 pb 


53if 5 pb 


502 


1.0±S;1 P b 


4.0±i;| pb 


H+3 Pb 


704 


0.10+° ;°1 pb 


0.56+°;?? pb 


2.o±8:g P b 


934 


10±| fb 


82±™ fb 


o.39tfS:iI pb 


1181 


i-o±8:| fb 


13±I fb 


84+^ fb 



TABLE IV: <7 e g for models I, VI from inclusive NLO QCD production. The error is the scale variation error. We use 
the NNLO PDF's of MSTW (see Section IV for details). 



A large amount of parameter space remains in each model that can be probed at LHC. The allowed 
regions of the model parameters permit large ~ pb production cross sections, and thus significant early 
LHC event rates that are not ruled out by current Tevatron searches or flavour constraints. This is illustrated 
in Fig. 4 for y/s = 7 TeV with the errors for the y/s = 7 TeV contours and values of the effective cross 
section for y/s = 10, 14 TeV given in Table IV. 



a e ff contours model I 



a e jf contours model VI 




> 

CO 
P. 600 




m 3 [GeV] 



m 3 [GeV] 



FIG. 4: Contours of constant rif<j determined from analytic NLO QCD pair production in models I, VI as described 
in the text. The labels on the plot correspond to y/s = 7 TeV. See Table IV for the errors and values for y/s = 
7, 10, 14 TeV for the shown contours of constant <j B g. In model I the ej mass correction is negligible for all three 
generations. We include in the contour the £3 correction to the vector-like quark mass dependence for model VI. For 
the first two generation vector-like quarks in model VI the contours are coincident with the left edge of the plotted 
contour and straight lines as the ei 2 mass correction is negligible. 
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V. CONCLUSIONS 

We have examined all SU(2)l singlet vector-like quark models with Q = —1/3 or Q = 2/3 that satisfy 
MFV while the spinor components are restricted to transform as a triplet under one of the flavour groups 
of Gf- We have directly solved these models and determined the mass spectra and mixing angles with 
the SM quarks. This has identified two of these models as naturally phenomenologically viable due to the 
predicted mixing angles of the model and the effective SM CKM matrix structure. These models are more 
phenomenologically viable than the remaining models due to their predicted mixing which depends on their 
group structure and representation content. However, we note that our analysis is done under the assumption 
that a perturbative analysis is appropriate and relies on the mixing that these models experience with the 
SM quarks at leading order in the MFV expansion. If rj > 1 so that flavour breaking insertions of the top 
Yukawa are not suppressed then a further analysis is required to determine the constraints on the mass scale 
of these models and the LHC discovery potential. 

The constraints on the two models were explored in detail using EWPD, Tevatron constraints and flavour 
observables and much parameter space remains for LHC to explore. We have demonstrated this with the 
effective NLO QCD production cross sections for the viable parameter space. 

Although the matter content we have studied has not been introduced to directly solve the hierarchy 
problem, such matter content may be a component of a new physics sector that does solve the hierarchy 
problem. The earliest evidence of such a sector at LHC could emerge from electroweak scale coloured 
states that satisfy known flavour constraints on the weak scale. 

We emphasize that if vector-like quarks are discovered in the early LHC era, they will likely have 
my < 1 TeV due to event rate limitations. In such a scenario the compatibility of the vector-like quarks 
with the constraints we have explored will be an important consistency check on a vector-like quark model 
consistent with an LHC signature. Due to their flavour structure, including the constraints of MFV, the 
models we have explored and constrained in detail are particularly promising SU(2)l singlet vector-like 
quark models for LHC discovery. 
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Appendix A: Remaining MFV model spectra and mixing angles 
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where 



Model II 



K f = gdi ni = 1,4 = 1, U{V, L) = U(V, R) = U(d, L) 



U 




w 



D 2 



(Mf y y o 

K) 2 + (mi) 2 + (Mf ys ) 2 




2g?g|+£f+2>/2&ftfe 
^+2g?gg-2>/2giftfe 






1 1,2,3 



M 



Phys_ ., kl + ZZUi-^^biz 



M 



e 2 +e 3 



(Al) 



In this model we identify the three eigenvalues directly proportional to M. d to be identical with the SM 
quark masses. The corresponding decay widths are obtained through rescaling the couplings of Eqs. 
(18)-(20). In this model, the mixing with the SM quarks is not naturally small unless one is in the 
decoupling limit rri2 — > oo. For this one and the following models in the Appendix, the results for the 



vector-like up quark models are trivial to obtain following the procedure discussed in Section II D 



Model III 



0, 4 = g d , 4 = g d , U(V, L) = U(d, L) , K(V, R) = U{d, R) 



( 



C+A 



U 



4653 



■c. 



4€ 2 ?3, 
C+- 



+ 1 



\ 



.4£ 2 5 3 , 

1 



+1 



\ 4653 



+1 



/ 



c_ 



1 



w 



c 



2^53 

1 



D 2 



*+* j \ w&) 2+i \/(^) 2+i / 

o {M' d f(e 2 +zi + \s+\ 



, 4 «2?3 J V V 4 C253 

r0\ 2 ( C 2 



where 



5 = J 4$ + (8£f - 4)£ 2 + (2^1 + I) 2 , C ±± = -2H ± 1H ± 6 + 1 . 



(A2) 



We identify three of these masses with the SM down quark masses and three of them with the new 



vector-like down quark states. We discuss the constraints on this model in Section HID 
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Model IV 



l,«2 = <?i,«3 = l. U(V,L)=U(d,R),U(V,R)=U(d,L) 



D 2 



mi) 2 + 0[{M 



0^21 







Mi 




+ 0[(M d f 



We refrain from including the mixing angles for this model in the Appendix due to their length and compli- 
cation, although we have determined them using the described procedure. The key point is that they do not 
provide naturally small non SM contributions to quark couplings to the W and Z. Also, as noted earlier, 
this model predicts non-hierarchical down quark masses and can only be viable as long as one chooses 
mf < rrid, where i = 1 or 3. 



Model V 



K l — Qdi K-2 ~ Qui K 3 — Qu 



D" 



U - 1 


1 £2 -6 




V> 6 



w 



£7. L>i 



1 



n 9 r< 2 1 .£ 2 u A/> phySN i 2 1 fA/> phySN i 2 €3+ 2 ^^+ 2 v / 2gi6g3 



^+^^-2^6 66 



where 



/"< _ g £2 ^3+y /2 €i £ 

° 4 - >fi 2 6 ef+«| 



2 C 

?1.2,3 



- 1 e 2 +e 3 



, M) 



(A3) 



In this model it is not possible to choose U(V, L/R) in such a way that the 3x3 subblocks of Eq. dT]) 
are diagonal. Nevertheless, we take U(V,L) = U(d,L), U(V,R) = U(u,R) and write the two left 
submatrices in a diagonal form and the two right submatrices as diagonal matrices multiplied by Vckm- 
We use the Wolfenstein parametrization and neglect all terms of order 0{\ 2 ) and higher. We assume that 
m d , m° are approximately equal to the SM physical masses up to small corrections, expanding in ra- as 
usual to simplify the eigenvalues and eigenvectors. The rotation matrices in this model are identical in 
form as in model II with the replacement S — > Q. The physics of this model is substantially the same, and 
unnatural with respect to flavour constraints. The only difference is that the light up quark masses receive a 



1.2 



large correction to their Higgs coupling as d 3> 1 due to M.u <C v. 
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